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The most general stationary black-hole solution of Einstein-Maxwell theory in vacuum is the Kerr- 
Newman metric, specified by three parameters: mass Af , spin J and charge Q. Within classical 
general relativity, one of the most important and challenging open problems in black-hole perturba- 
tion theory is the study of gravitational and electromagnetic fields in the Kerr-Newman geometry, 
because of the indissoluble coupling of the perturbation functions. Here we circumvent this long- 
standing problem by working in the slow-rotation limit. We compute the quasinormal modes up to 
linear order in J for any value of Q and provide the first, fully-consistent stability analysis of the 
Kerr-Newman metric. For scalar perturbations the quasinormal modes can be computed exactly, 
and we demonstrate that the method is accurate within 3% for spins J/Jmax J; 0.5, where Jmax 
is the maximum allowed spin for any value of Q. Quite remarkably, we find numerical evidence 
that the axial and polar sectors of the gravito-electromagnetic perturbations are isospectral to linear 
order in the spin. The extension of our results to nonasymptotically flat space-times could be useful 
in the context of gauge/gravity dualities and string theory. 

PACS numbers: 04.70.Bw, 04.25. Nx, 04.30.Db 



Introduction. In Einstein-Maxwell theory, black holes 
(BHs) that are stationary, asymptotically flat end-states 
of gravitational collapse must be axisymmctric [l| . Clas- 
sic uniqueness theorems reviewed in [^ show that regular, 
stationary electrovacuum BH space-times in four dimen- 
sions are described by the Kerr-Newman (KN) metric [j] , 
characterized by mass M, angular momentum J and elec- 
tromagnetic charge Q. When Q = the KN solution 
reduces to the Kerr metric, and for J = it reduces to 
the Reissner-Nordstrom (RN) metric. When both Q and 
J are nonvanishing the space-time is endowed with an 
induced magnetic field, and its magnetic dipole moment 
corresponds to the same gyromagnetic ratio g = 2 as the 
electron Q. This observation led to some speculation 
that the KN metric could be used as a classical model 
for elementary particles (see e.g. Q). 

Charge is unlikely to play a significant role in astro- 
physics [a, 13 : but the KN metric is still a precious theo- 
retical laboratory to investigate Einstein-Maxwell theory 
in curved space-time. For this reason the linearized dy- 
namics of test fields on a KN background have been in- 
tensively studied in the past. The scalar Q, neutrino Q, 
massive spin-1/2 [i3,[ll| and Rarita-Schwinger [T^l equa- 
tions in the KN metric can all be solved by separation of 
variables. The scattering of charged scalars and fcrmions 
in near-extremal KN space-times recently acquired spe- 
cial interest in the context of the KN/Conformal Field 
Theory (CFT) conjecture [Hill. 



The KN space-time is one of the simplest prototypes of 
the interplay between matter and curvature summarized 
by Wheeler's famous statement that "matter tells space- 
time how to curve, and space-time tells matter how to 
move." Despite their importance, theoretical investiga- 
tions of the interplay between gravitational and electro- 
magnetic perturbations in the KN metric are still in their 
infancy. The reason is a major technical stumbling block: 
most methods to compute quasinormal modes (QNMs, 
see |15l - [l8j for reviews), greybody factors and scatter- 
ing amplitudes rely on separability, and despite several 
attempts |19l - [2l| , at present no one has been able to sep- 
arate the angular and radial dependence of the gravito- 
electromagnetic eigenfunctiontl^. 

This work is the first consistent analysis of gravito- 
electromagnetic perturbations of KN BHs. We circum- 
vent the decades-old coupling problem using a recent 
framework to study generic perturbations of spinning 
BHs in the slow- rotation limit (23, 01, which is based 
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^ The last chapter of Chandrasekhar's monumental 1983 mono- 
graph [2a is dedicated to an incomplete treatment of this prob- 
lem. Quoting from |22| : "It does not appear that the methods 
developed [...] for the treatment of the gravitational perturba- 
tions of the Kerr black-hole can be extended in any natural way 
to the treatment of the coupled electromagnetic and gravita- 
tional perturbations of the KN black-hole. The origins of this 
apparently essential difference in the perturbed Kerr and KN 
space-time may lie deep in the indissoluble coupling of the spin- 
1 and spin-2 fields in the perturbed KN space-time - a coupling 
which it was possible to break only for very special reasons in 
the perturbed RN space-time." 



on a similar approach used in the past to study slowly 
rotating compact stars |25l - [28| . We summarize here some 
of our most interesting findings: (i) We present the first 
self-consistent calculation of scalar, electromagnetic and 
gravitational QNMs of the KN metric, (ii) Since none 
of these modes is unstable, our calculation provides solid 
evidence for the stability of the (noncxtrcmal) KN met- 
ric, (iii) In the scalar case we can compare our results 
to an exact calculation, that does not rely on the slow- 
rotation limit. We find that the perturbative analysis is 
valid when J/Jmax ^ 1 (where Jmax is the maximum 
allowed spin for any given Q) and that scalar QNM fre- 
quencies are accurate within 3% for spins J/ Jmax ^ 0.5, 
which suggests a similar level of accuracy for the gravito- 
electromagnetic modes, (iv) Last but not least, we find 
the remarkable result that axial and polar QNMs (corre- 
sponding to perturbations that have odd or even parity, 
respectively) are isospectral to linear order in the spin. 
Formalism. The KN metric is the most general sta- 
tionary electrovacuum solution of Einstein-Maxwell the- 
ory. Its full form in Boyer-Lindquist coordinates can be 
found, e.g., in [22|. Here and in the following we linearize 
all quantities in the spin parameter 5, = a/M = J/M'^ 
(in geometrical units G = c = 1), neglecting terms of 
order O(a^). To this order, the KN metric reads 

dsl ^ -Fdt^ + F-^dr^ - 2msm'^ Mipdt + r^d'^fl , (1) 



where F(r) = {r-r^)ir-r+)/r^, r± = M±yjAP - Q^ 
arc the horizons of a RN BH. the gyromagnetic term is 

Tu{r) = 2hM'^/r - dQ'^M/r'^, (2) 

and the background electromagnetic potential is given by 
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Note that the presence of both rotation and charge {dQ ^ 
0) induces a magnetic field in the (?9, Lp) directions. 

We derive the equations describing gravito- 
electromagnetic oscillations in the slow-rotation 
approximation by linearizing the Einstein-Maxwell 
equations with respect to both the oscillation amplitude 
and the BH spin parameter a, and by expanding the 
perturbations in a complete basis of tensor spherical 
harmonics. As a consequence of using this basis in 
a nonspherical background, the linearized equations 
display mixing between perturbations with different har- 
monic indices and opposite parity J23l - l25l l27t . However, 
the latter do not contribute to the QNM spectrum to 
first order in 5 [13, [H, [H . 

Our main analytical result consists of two sets of cou- 
pled, second-order equations (one for the axial and one 
for the polar sector, respectively) which fully describe 
gravito-electromagnetic oscillations of a KN BH to first 
order in the spin. In the frequency-domain, and assuming 
a time dependence e^*'^*, they read (schematically) 
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where i,j = 1,2, i ^ j (there is no sum over the 
indices i,j), we have defined the differential operator 
V = dl^ + uj^ - F£{i+ l)/r^ and r* is the standard 
tortoise coordinate, such that dr,r ~ F. The functions 
Z~ and Z^ are linear combinations of axial and polar 
variables, respectively, and they are also combinations of 
gravitational and electromagnetic perturbations. 

The explicit form of the axial and polar potentials 
]/(^±) and M^(*'^) is quite formidable. It will be pre- 
sented in the accompanying paper [30 1 and in a publicly 
available Mathematica notebook [3l[ . What matters is 
that Eqs. (|3]) display the same symmetries as the master 
equations for a RN BH jg^] , and indeed they exactly re- 
duce to the latter in the nonrotating case. In addition, 
Eqs. (m contain two first-order corrections in a. The 
first term is responsible for a Zeeman-like splitting of the 
eigenfrequencies, which breaks the degeneracy in the az- 
imuthal index m. The second line in Eq. ^ is more 
interesting: this term couples the function Z^ with the 
function Z2 , and the function Z^ with the function Z^ . 

Once physically-motivated boundary conditions are 
imposed. Eq. ^ defines an eigenvalue problem for the 
complex frequency uj = ujr + iuji. The boundary condi- 
tions for QNMs read simply 
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The near-horizon solution displays the typical frame- 
dragging effect occurring near spinning BHs, where 
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is the angular velocity at the horizon of locally nonrotat- 



ing observers, and an 



Jmax/M2 = ^1 - {Q/Aiy 



is 



the maximum spin parameter of a KN BH 
Numerical Results. The numerical solution of the ax- 
ial and polar perturbation equations ^ is challenging, 
because their explicit form is very complicated 30|, |31 



Robust numerical methods to solve the coupled eigen- 
value problem given by Eqs. Q with the boundary condi- 
tions (O are reviewed in [S^] . We have integrated the cou- 
pled system Q and computed the corresponding eigen- 
frequencies using two independent methods: a highly- 
efficient matrix- valued continued fraction technique and 
direct integration [SG . |33 | . When both methods are ap- 
plicable they validate each other, in the sense that the 
results agree within numerical accuracy. 

For any given Q, our analysis allows us to extract the 
first-order corrections to the complex QNM frequencies: 
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where ui^ j are functions of Q and of the multipolar index 
E, whereas the m-dependence has been factored out [2J|. 
As a test of the slow-rotation approximation, we have 
computed the scalar QNMs of a KN BH to first order in a. 
These modes can be computed exactly in the Teukolsky 
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FIG. 1. Top panel: percentage error A = 10^|1 — Wsiow/ttJcxactl 
induced by the slow-rotation approximation in the real part 
of the fundamental £ — m = 1 scalar mode. Bottom panel: 
percentage error for the imaginary part of the same mode. 
The errors are only mildly sensitive to Q if plotted as functions 
of J/ Jmax, where Jmax is defined below Eq. © . Similar results 
hold for other scalar modes [3011. 



formalism JSJI, so they give us the precious opportunity 
to estimate the errors introduced by the slow-rotation 
approximation. For any stationary and axisymmetric 
space-time, the scalar modes at first order in the angu- 
lar momentum arc governed by a master equation 24 1 
whose corresponding eigenvalue problem can be solved 
with standard continued-fraction techniques |34l . |35| . 

In Fig.[T]wc show the relative error of the slow-rotation 
approximation with respect to the "exact" result, com- 
puted by solving the scalar equation in a KN background 
via continued fractions [33|. In particular, the top (bot- 
tom) panels show the percentage deviation A for the 
real (imaginary) part of the fundamental £ = m = 1 
scalar mode at fixed values of Q. The slow-rotation ap- 
proximation is accurate within one percent as long as 
J/Jmux ^ 0.3, and it is still accurate within 3% for 
J/ Jmax S 0-5- Similar results also hold for other val- 
ues of £ and m, and for the first few overtones [30|- Note 
the near-universal behavior of the percentage errors as 
functions of J/ Jmax = a/omax for all values of Q. In- 
deed the parameter flmax, which appears explicitly in the 
QNM boundary conditions Q , plays a fundamental role 
in our pcrturbativc scheme and the slow-rotation approx- 
imation is accurate only far from cxtremality, i.e. when 
Vl - QVM2 < 1. 



a < Omi 

Figure [5] shows our main numerical results for the 
fundamental gravito-electromagnetic perturbations with 
£ = 2, the most relevant for gravitational-wave emission 
(see e.g. [36|). In each panel we show four curves, corre- 
sponding to the axial and polar "gravitational" and "elec- 
tromagnetic" modes (as defined in the decoupled limit, 
Q — > 0). The zeroth-order terms shown in the left panels 
are simply RN QNMs; they agree with continued-fraction 
solutions [3J] of the equations first derived by Zcrilli [33] ■ 
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FIG. 2. Zeroth-order (left panels) and first-order (right pan- 
els) terms of the slow-rotation expansion of the KN QNM 
frequencies [cf. Eq. ((Tjl]. All quantities are plotted as a func- 
tion of Q/M , and they refer to the fundamental mode (n = 0) 
with £ = 2. The lower part of each panel shows the percent- 
age difference between axial and polar quantities: our results 
are consistent with isospectrality to 0(0.1%) for the real part 
and to 0(1%) for the imaginary part of these modes. 



We carried out a more extensive QNM calculation 
working in the axial case, where our results can be ver- 
ified using two independent methods. By virtue of the 
isospectrality between axial and polar modes visible in 
Fig. [2] and discussed below, these results cover the whole 
QNM spectrum of slowly rotating KN BHs. We found 
that the zeroth- and first-order quantities shown in Fig. [2] 
(plus analogous calculations for £ > 2 and the first over- 
tones [3^ ) are well fitted by functions of the form 
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Here we have defined a parameter y — 1 — Omax, which 
is in one-to-one correspondence with Q/M , but is better 
suited for fitting. The coefficients fi for the fundamental 
£ = 2 gravito-electromagnetic modes are listed in TableU) 

Stability. None of our numerical searches (for < Q < 



M, J < J„ 



and 



2,3,4) returned exponentially 



growing modes. This confirms early arguments by Mash- 
hoon in favor of the stability of the KN metric [S^ . Mash- 
hoon's results apply only to the cikonal limit (.^ ^ 1) 
and they rely on a somewhat heuristic geodesic analogy, 
rather than on a self-consistent treatment of the pertur- 
bation equations. In this sense, our QNM calculations 
provide the first, fully consistent numerical evidence for 
the stability of the KN space-time. 

Isospectrality. Gravito-electromagnetic perturbations 
of Schwarzschild and RN BHs in general relativity have a 
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/o /l /2 /a /4 


4°' (2,0,1) 
J^'> (2,0,1) 
^T (2,0,1) 
4" (2,0,1) 


0.4576 0.2659 0.0118 0.1228 -0.1382 
0.0712 0.0769 0.0596 0.0727 -0.0216 
-0.0950 -0.0184 0.0137 0.0132 0.0107 
0.0007 0.0043 0.0060 -0.0089 0.0366 


u^°^ (2,0,2) 
^k'^ (2,0,2) 
4"^ (2,0,2) 
4" (2,0,2) 


0.3737 0.0525 0.0607 -0.0463 -0.0070 
0.0628 0.0676 0.0209 0.0823 -0.0810 
-0.0890 -0.0055 0.0024 0.0214 -0.0084 
0.0010 0.0014 0.0091 0.0174 0.0145 



TABLE I. CoefHcients of the fit ^ for the real and imagi- 
nary part of the fundamental (n = 0) gravito-electromagnetic 
modes with i = 2. We denote by s = 1 and s = 2 the 
modes that in the decoupled limit, Q — ,- 0, are electromag- 
netic and gravitational, respectively. The fits ([8]) reproduce 
the data to within 1% for uj] , and to within 0.1% for the 
other quantities, for any Q < 0.95M. Similar fits have com- 
parable accuracy also for (. > 2 and for the first overtone [301 . 



noteworthy property, first proved by Chandrasekhar [22l | : 
even though the polar and axial sectors of the pertur- 
bations are described by completely different potentials, 
their QNM spectra are identical [l3]. Mathematically, 
this happens because the polar and axial potentials can 
be written in terms of a single "superpotential" (cf. Sec- 
tions 26 and 43 of [221 ) . This property can be interpreted 
as being due to "supersymmetry," in the sense of nonrel- 
ativistic quantum mechanics |39l44l| . 

Isospectrality is easily broken: e.g., it does not hold if 
the cosmological constant is nonzero 143-1441 . if the under- 
lying theory is not general relativity [451 146| , or in higher 
dimensions (cf. Appendix A of [131 )■ The left panels of 
Fig. [2] (which refer to the RN limit) show that polar and 
axial modes are isospectral within our numerical accu- 
racy. Given the complex form of the polar equations, 
this is a nontrivial check of our numerical techniques. 

A priori, there is no reason why such a remarkable 
and fragile symmetry should hold true also for rotating 
(KN) BHs. A tantalizing result of our numerical study 
is strong evidence that the axial and polar sectors of KN 
gravito-clcctromagnctic perturbations are indeed isospec- 
tral to first order in the BH spin. The left panels of Fig. [2] 
show that the linear corrections w]j \ are identical func- 
tions of Q for axial and polar modes within the numerical 
errors (which are dominated by uncertainties in the di- 
rect integration used to compute polar modes). Various 
arguments can be made to support the claim that the 
observed deviations from isospectrality arc of a purely 
numerical nature: 1) isospectrality is verified to a higher 
level of accuracy far from extremality: this is consistent 
with the fact that QNMs are more challenging to com- 
pute in the extremal limit; 2) the deviations from isospec- 
trality shown in Fig.[2]are roughly constant or decreasing 
functions of Q (at least for Q < 0.8A/) and they are af- 
fected by a small residual error even when Q — Q, where 
isospectrality must hold exactly; 3) the direct integration 



method is more accurate as I grows and, correspondingly, 
the deviations between axial and polar modes decrease; 
4) finally, we verified that the error can be reduced by 
increasing the accuracy of the integrator. 
Outlook. It is tempting to conjecture that the isospec- 
trality we found at linear order may in fact hold exactly, 
at all orders in rotation. In order to verify this hypoth- 
esis it will be crucial to include effects of second order 
in the spin - a formidable undertaking. At second order 
the causal structure of a spinning metric starts differing 
from the nonspinning case, and parity-mixing terms ap- 
pear in the perturbation equations |24{ . If isospectrality 
were to hold true also at second order, there would be no 
fundamental reason to believe that it should be broken 
at higher orders. However, let us stress that isospectral- 
ity is a highly nontrivial property even at linear order 
in rotation, in view of the mixing of gravitational and 
electromagnetic perturbations. Hopefully our work will 
stimulate further studies to verify whether isospectrality 
is an exact property of the KN space-time. Besides brute- 
force extensions of our work to higher orders in rotation, 
other possible means of study ing this problem include nu- 
merical time evolutions (cf. |47H49f) or analytical work, 
perhaps along the lines of [4l| . 

Other interesting extensions of this work concern 
asymptotically anti-de Sitter (AdS) space-times. Our ap- 
proach can be easily applied to compute the QNMs of 
(slowly rotating) KN-AdS BHs. This would be useful in 
the context of the AdS/CFT correspondence [501, "which 
predicts that these solutions are dual to thermal states 
of a CFT living in a rotating Einstein universe [5l|, [52l . 

Such an extension would also be interesting in the con- 
text of supergravity. To clarify this point, let us recall 
that the QNMs of asymptotically flat, extremal RN BHs 
have a remarkable property: electromagnetic perturba- 
tions with angular index (. are isospectral with gravita- 
tional perturbations with index (. -\- 1 [5j1. This is re- 
lated to the fact that, when embedded in A^ = 2, four- 
dimensional supergravity, these solutions preserve part 
of the supersymmetry. Using this property, it is possible 
to prove that the one-loop corrections to the BH entropy 
cancel [SJl. In the case of rotating, asymptotically fiat 
BHs this reasoning does not apply, since these solutions 
are not supersymmetric. However, certain KN-AdS BHs 
embedded in A'^ = 2 four-dimensional supergravity, pre- 
serve half of the supersymmetry [55|,|56| . Since these solu- 
tions can be slowly rotating, the techniques developed in 
this paper could be used to extend the arguments of [5J1 
to supersymmetric KN-AdS BHs. 

Last but not least, it would be interesting to extend our 
calculation of KN QNMs in the context of the KN/CFT 
conjecture [13,113, "which predicts that the QNMs of the 
near- horizon KN geometry correspond to the poles of the 
retarded Green's function of the dual chiral CFT [531 • 
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